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400 399 398 397 396 395 394 393 392 391 590 389 388 387 386 385 3584 383 382 381 
525 324 323 322 321 320 319 318 317 316 315 314 313 312 311 310 309 308 307 380 
326 257 256 255 254 253 252 251 250 249 248 247 246 245 24H 243 242 241 306 379 
327 258 197 196 195 194% 193 192 191 190 189 188 187 186 185 184 183 240 305 378 
328 259 198 145 144 143 142 141 140 139 138 137 136 135 134 133 182 239 304 377 
329 260 199 146 101 100 99 98 97 96 95 94% 93 92 91 132 181 238 303 *76 
= 330 261 200 147 102 65 64 63 62 61 60 59 58 57 90 131 180 237 302 375 
331 262 201 148 103 66 37 36 35 34 33 32 31 56 89 130 179 236 301 374 
332 263 202 149 104 67 38 17 16 15 W 13 30 55 88 129 178 235 300 373 
333 264 203 150 105 68 39 18 5 4 3 12 29 54 87 128 177 234 299 372 
334 265 204 151 106 69 40 19 6 1 2 11 28 53 86 127 176 233 298 371 
335 266 205 152 107 70 41 20 7 8 9 10 27 52 85 126 175 232 297 370 
336 267 206 153 108 71 42 21 22 23 2h 25 26 51 84 125 17% 231 296 369 
337 268 207 154 109 72 43 4Y 45 46 47 4B 49 50 83 124 173 230 295 368 | 
338 269 208 155 110 73 74 75 76 77 78 79 80 81 82 123 172 229 294 367 
339 270 209 156 111 112 113 114 115 116 117 118 119 120 121 122 171 228 293 366 
340 271 210 157 158 159 160 161 162 163 164 165 166 167 168 169 170 227 292 365 
341 272 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 291 364 
342 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 363 
@ 343 34H 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 
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20 x 20 


On the front cover of this issue is a 20 x 20 array 
of the numbers from 1 to 400 arranged in a spiral. The 
cells in the array are identified by columns, from one to 
20, left to right, and by rows, from one to 20, top to 
bottom. Each cell can thus be labelled as alr,J : Thus 
we wish to have: 


A(1,1) = 400 
A(2,2) = 324 : 
A(3,3 = 256 ee 
° x a 
a e = 
A( 18,18) = 236 d 
A(19,19) = 290 @ 
A(20,20) = 362 . 


Part 1 of this Problem, then, is to generate the 
400 numbers in the proper cells in the array A. 


We then wish to be able to shift any row or column 
by any number of cells, in circular fashion. That is, we 
want a subroutine to shift row Q, R cells to the right, with 
the numbers shifted off the array at the right appearing in [ 
the same row on the left end. Thus, for Q = 11, R = 3, 
row 11 will become: 


eaceemeos sf 9334 265 ........ 53 86 le7 42176 
If R = 20, there will be much activity, but the row will 
appear unchanged. To shift a row to the left by W places, 


let R = (20 - W). 


Similarly, we want a subroutine that will shift 
column S up by T cells. Again, if T = 20, the column 


will appear unchanged. To shift a column of the array down 
by Z, let T= (20 - Z). 
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The center of the array is this: 


The following actions are to be performed: 


Shift row 11 1 place to the right 
Shift column 11 1 place up 

Shift row 10 1 place to the left 
Shift column 10 1 place down. 


Following that, the rows and columns just outside the 
center (namely, rows 9 and 12; columns 9 and 12) are to be 
similarly shifted two places; that is, 


Shift row 12 2 places to the right 
Shift column 12 2 places up 

Shift row 9 2 places to the left 
Shift column 9 2 places down. 


This procedure continues, moving out from the center, until 
finally we: 


Shift row 20 10 places to the right 
Shift column 20 10 places up 

Shift row i 10 places to the left 
Shift column 1 10 places down. 


The final status of the array is shown on the back cover. 


Producing that result is Part 2 of this Problem. 


Everything called for is now done, and the result is 
recorded. Precisely; and we now have a splendid coding 
assignment, readily defined, with an enormous amount of 
detailed work to be performed--all leading to a known 
result (that is, it 1s known if we have done our work 
correctly). fry it--it appears to be quite simple. 


ae) 
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Take/Skip Revisited 2 


Probably the most spectacular problem we have so far 
published (and certainly the one that attracted the most 
interest and work) first appeared as the K-Level Sieve in 
issue 38, but was soon dubbed the Take/Skip Problem. A 
brilliant solution by two University of Toronto students, 
Tom Duff and Hugh Redelmeier, appeared in issue 43. 


Contributing Editor Edward Ryan observed that the 
technique used in "...Or Not Recurse" in issue 75 could 
apply to the Take/Skip problem, and that therefore a fairly 
efficient implementation of Duff and Redelmeier's solution 
could be made in elementary BASIC, Let us restate the 
original problem: 


Start with the positive integers. At each level, K, 
of the procedure to be followed, the numbers that 
are outputted from level K-l are to be sieved by 
accepting K numbers and rejecting K numbers, 
alternately. What numbers will survive all 

levels of sieving? 


The positive integers form level zero. Level 1 
accepts one number and rejects the next; thus the 
output of level 1 (which constitutes the input to 
level 2) is the sequence of odd numbers. 


Subsequent levels will deal with the following sequences: 

2 1 3 9 il 17 419 5) eT 33 35 Al 43 
Chmesyten 33) 40% obi 2 Si aieas oe Commoner 
Oyeeoer oi 73 gT5 gS) 123 mae eG 
9 


3 1 
4 1 
5 1 25 57 145 147 193 #195 201 ... 


3 
3 
3 


Each level passes on to the next level exactly half 


of the numbers it receives. Thus, for a million numbers 
treated at level zero, just one number will emerge at 
level 19. 


The point to focus on, however, is that each level 
treats just one number at a time, and rejects it (and that 
number is then seen no more) or accepts it and passes it 
on to the next level. So if we observe the numbers that 
survive to any given level, say level L, then as long as 
those numbers are less than 2 to the power L, they are 
numbers that will survive all levels of the sieve. 


G--9Ldd 


apebiig yayong = |1S/DVL 


*TSAaT 48U3 9B e9Bp aug BuyQe1aues ‘Z TaAaT 
4B BupTyJtom sft *Z voUuadayTer ye Buy qaegs ‘auTgnod sTuy 


S[NS8I 3yQ 2aS 09 USTM NOA TaAaT BUM 4B BaUTJap % 
03 08 04 sSTaAaT JO aAsqumu 94, sauTJep J 

*TAAST 09 [eAeT WOIT KX BuTAow soz aeHONG @ BT OO 
[TOIZUOD YOATMS & ST W 

B4ep 949 yo aJoYeaauad ayyq BST xX 

ae | ty T9AaT 4B ‘Bets e& ST @ yORA 


WY Teast BZuypucdseiiz09 ayy Jos Jequnoo eB s}t gq yory 
aequnu [T@eAeT aya ST ¥ 


xeu) suTgnod 
Terauas 
3syuy 903 9H 


04 T=HI dog 
T= (1)a 
(6) (I)a 39S 
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General routine 
for level K 
starts here 


):0>= 
Jee 
# 


“F(K 


B(K)+1—> B(K) 


So all we need is a bunch of counters and flags. 


The flowcharts given here have been made to operate through 


level 50 (and 1125899906842624 is the 50th power of 2). 


We set up a counter, B, for each level, and all the 
B's are initialized to zero. We also establish a flag, 


at each level, all initialized to one (standing for "accept') 


The mechanism at Reference 2 generates the level 2 output. 


FF 


Beginning with level 3 (K = 3), all levels operate the same 


way, as shown in the second flowchart. 


Consider level 10. It receives numbers, one at a 
time, from level 9 (transmitted via word CC). Starting 
from the beginning, it should pass 10 numbers on to level 
11, then reject the next 10 numbers it receives, then pass 
on the next 10, and so on. Its counter, B(10), counts 
up to the level number 10, and then resets to zero to get 
ready to count to 10 again. Its flag, F(10), flipflops 
from zero to one to zero to indicate whether it is in the 
accept (1) or reject (0) mode. When any level is in 
reject mode, the contents of CC is made zero, and control 
returns to Reference 2, to generate the next number at 
level 2. 


Even in interpretive BASIC, this procedure is 
reasonably fast. A possible program is included. 


kKeekkkkkkKkkKkeenenen 


Most Back Issues Are Still Available: 


eat. 


APR 
DEC 


al = B 
zaos 8 
Ieee 237 cal 45 x 
mK 19 
Sls cesS Vol. 


8 

KG gs Noten 
HOM dia Ie Nf 1h 3 21 Vol. 
ome 24 625) 26 «27» 28° 229 
34 9{ 36 37 38 39 4O Al 
46 47 48 49 50 51 52 53 
58 59 60 61 62 63 64 65 


Com alae Tae (4 675 Vol. 


43° 44 45 Vol. 
55 /50m54) Vol. 
67 68 69 Vol. 


KR 

(= 

~] 
pee oR pK SEF 


Nn Ww F WwW PD 


« 


1973 
1974 
1975 
1976 
1977 
1978 
1979 


PC76--7 


PC76--8 


100 DIM B(50) 

110 DIM F(50) 

120 FOR I = 1 T 50 

130 ae = QO 

140 F(I) = 1 

150 NEXT I 

160 X = -5 

170 M = 0 

500 IF M = O THEN 600 

5lo X=X+2 

520 M = O 

530 G@ TH 700 

600 M= 1 

610 X= X+t+6 

70O CC = X 

710 K = 3 

1000 IF F(K) = O THEN 1100 
1010 IF K > 10 THEN 1200 
1020 B(K) = B(K) + 1 

1030 IF B(K)>=K THEN 1300 
1040 B(K) = 

1050 IF F(K) = O THEN 1400 
1060 anes = 0 

1070 G% TS 1300 

1100 cc = 

1110 B(K oe = B(K) + 1 

1120 IF B(K) < K THEN 500 
1130 B(K) = 0 

1140 IF F(K) = O THEN 1500 
1150 F(K) = 0 

1160 @% T% 500 

1200 IF CC = O THEN 1020 
1210 PRINT K, CC 

1220 GG T% 1020 

1300 K =K+1 

1310 IF K > 48 THEN 500 
1320 GS T@ 1000 

1400 F(K) = 1 

1410 @g T% 1300 

1500 F(K) = 1 

1510 @% TS 500 


Possible program in BASIC for the TAKE/SKIP Bucket Brigade 


Line 500 is Reference 2 of the flowchart 


Line 700 is Reference 3 of the flowchart 
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Penny Flipping VII 


There are two stacks of coins numbering C coins in 


each stack. Initially, both stacks are all heads up. 
The accompanying illustration is for C = 3. 


The first stack follows the procedure of the first 
Penny Flipping problem; namely: 


Flip the top coin, then the top two, then 
the top three,..., and so on until the 
entire stack is flipped; then start over 
with the top one, the top two,... 


The number of coins that are tails up in the first stack at 
each stage is the number to be flipped in the second stack. 
Thus, the top coin of the first stack is flipped to start, 
which gives one tail, which dictates flipping one coin of 


the second stack. This procedure continues until both 
stack again become all heads up. For C = 3, this takes 
36 flips, as shown. The results for cases 2 through 40 
are given. 


The results for previous penny flipping problems 
were distributions that oscillated wildly. This latest 
version, due to Ralph Montgomery of St. John's Community 
College, which would appear to be even wilder from its 
description, turns out to be (at least for 40 cases) 
remarkably well behaved. 


As usual, the point of the proodlem is to provide an 
interesting coding exercise (the penny flipping problems 
lend themselves well to integer BASIC or machine language) 
for which test results are available. We solicit further 
results for values of C greater than 40. 


700101210111200101011100132010121001010-: 
11100120001110010001110010001110010 
0001110000001110000001110000001110 
101012110120101011011011000100100100 
1110011110001100000000011113211211100 
oo000e0O0000000017212121212112112113212113a3ia 00 
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8078400 


Penny flipping VII results March 1979 
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2 THROWBACK Revisited 


The THROWBACK problem (Number 193) was on the cover 
of issue number 55. The problem was this: 


Given all the positive integers starting with 3: 


Som oa GameOee 10-37, ~ Veh. Ta eerie 
The leading 3 specifies throwing it (the 3) back 3 places, 


putting it between the 6 and 7. This process now 
repeats, as shown here: 


Seo 5}. mance) 10 11 12 13 14 
6 Sele /; k 8 9 10 11 12 13 14 
3 a 8 5 9 10 ii 12 13 14 


geet FORO ND .19- 6 LO; | 2) oS 


and we count the number of throws to bring each higher 
number to the leading position, giving a table: 


4 1 
5 2 
6 3 
it e) 
8 7 
9 10 
10 14 
19 


In issue 58, this table was extended to line 48, largely 

as a result of work done by William Bourn, who did it with 
& program in APL. The current interest in BASIC prompts 
us to present a program in that language to implement a 
solution. But this leads us to a new problem. 


In the original version, the results appear to follow 
the rule: 


Fo. = (4/3 ° Pea 1], 


Se ntl 
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at least for the first 48 results. This indicates that 
the number 100 will first appear at the head of the stream 
of numbers after some 2,731,288,000,000 throws have been 
made. 


Now suppose that the starting sequence of numbers 
began with 10, rather than 3. In the BASIC program, these 
two changes will get the result empirically: 

120 A(I) =I+9 
140 X = 11 


and the resulting table begins: 


21 au 19 9 27 26 
12 2 20 10 28 29 
13 3 21 12 29 32 
14 4 22 14 30 36 
15 5 23 16 31 40 
16 6 24 18 32 AS 
17 i 25 20 33 50 
18 8 26 23 34 56 


The new problem is to reach a logical conclusion 
as to when the number 100 will first arrive at the head 
of the stream. 


100 DIM A(100) 

110 FOR I = 1 TH 100 

120 A(I) =I+2 

130 NEXT I Program in BASIC for 


l4o x =4 
150 N=0 the original THROWBACK 


200 K = A(1) problem. 
210 FPR IT=1 WK 
220 A(I) = A(I+1) 
230 NEXT I 
240 A(Kt1) = K 
250 N=N+]1 
260 IF A(1) = X THEN 500 
270 GSTS 200 


500 PRINT X, N 
5lo X=X+t1 
520 GSTS 200 


530 END 
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Hopefully, a final report “any Is 

on the first Penny Flipping lip ’ 
Problem. See previous Ping 
reports in issues 23, 25, Ay ’ 
71, and 73. Wn 


The last reference (issue 73, page 18) reported on 
research by David E, Ferguson, he who runs the Ferguson Tool 
Company. Mr. Ferguson adds the following results now: 


2*( moa a) can be computed in at most 2k adds (average 3k/2 


adds). Since (2C+1) is always even, 


(2C+1)/2 
2 ( ) = +1(mod 2C+1) 


and since #(2C+1)<2C, n<C, n can be computed in at 
most 2C adds (average 3C/4 adds). 
Also, if 2Ctl is prime and 
(a) C=1 or C=2 (mod 4) f£(C) = C/d-1 where d is 
a proper odd divisor of C. 
Corollary 1: If C is a prime=1 (mod 4), then 
f(6) = Cc? -1 
Corollary 2: If C is twice a prime, f(C) = Com 
(b) C=3 (mod 4) f(c) = c°/d where d is a proper 
divisor of C. 
Corollary 3: If C is a prime= 3 (mod 4), then 
t(c) = cf 
(c) C=O (mod 4), either 


(c1) f(c) = 07 /a-1 where d@ is a proper 
even divisor of C; or 
(c2) f(C) = Cd where d is a proper odd 
divisor of C. p 
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(ad) f(C) 1s of the form nC-1 if and only if f(p,) is 
of the form nyP,-1 for every prime divisor, Py> 
of 2Ctl. (This accounts for the thinning out 
of this form for higher values of C.) 

(e) A general expression for f(C) can be written in 
terms of the prime factorization of 2C+1 which 
will work in all but very rare cases; namely, 
when 2P-+ 23 (mod p“); e.g., when p = 1093. 

(f) My "£(C)" 1s your "N"; my "n" in as +1 is not 


related to your "N", 


Meanwhile, using very little analytics, but clever 
coding (by Associate Editor David Babcock) and enormous 
amounts of computing power, the table of results for the 
first penny flipping problem from issue 71 can now be 
extended. For stacks of coins from 240 to 539, and from 
1000 to 1103, the number of flips is given in the accompany- 
ing table. 


In the review of the TRS-80 by Larry Clark in 


issue 75, Mr. Clark was credited with appearing 


in the film "Joss." The confusion was 
natural, inasmuch as he was in charge of the 
JOSS project at The RAND Corporation for some 
time. The film appearance, however, was in 


the AFIPS film "It's Your Move." 


2ho 


250 


260 


270 


280 


290 


300 


310 


320 


330 


340 


350 


360 


370 


380 
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390 


400 


410 


420 


430 


128188 


440 


450 


460 


470 


480 


490 


500 


510 


520 


530 


265188 


1000 


1010 


1020 


1030 


308000 
143142 
200400 
222666 
421680 
202004 
60360 
60420 
169343 
48432 
325220 
412488 
546480 
1026168 
1028195 
686140 
484631 
2G A000 
48864 
1038361 
79559 
346118 
104243 
11253 
11263 
897900 
1052675 
69836 
452320 
144060 
234840 
1062961 
359136 
161148 
1069155 
37260 
119139 
850340 
342540 
93510 


1040 


1050 


1060 


1070 


540799 
1083680 


1080 583199 
110262 
T7904 
1061340 
26016 
1080660 
141179 
152180 
505920 
395306 
1090 263780 
1139004 
432432 
796796 
1083060 
170820 
61376 
320324 
1113371 
268156 
1100 38500 
404066 
92568 
1216609 
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Palindromic Numbers 


As part of a larger problem, it is important to 
know whether or not the decimal representation of a 
number is a palindrome. A palindrome is something 
which "reads' the same forwards and backwards. 


For example, the numbers: 
1 44 323 5775 and -848 
are palindromes, while the numbers: 
12 344 67 -378 and 5050 
are not. 


The program segment below and the accompanying 
flowchart perform the required task. 


INTEGER PROCEDURE PALINDR@ME (N,B); 
INTEGER N,B; 


BEGIN INTEGER R,S,T; 
R:=0; 
S:=ABS(N); 
T:=2; 


D6 BEGIN 
R:= R*B + MOD(S,B); 
S:= S/B; 


IF R<>ABS(N) THEN PALINDROME := 1 
ELSE PALINDROME := T; 
END; 
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Number to be tested 


N= 
B = base 
ES T igs returned as: 
we Fe 1: Not a palindrome 
ost 2: Palindrome, no center 
DEES marker 
vog 3: Palindrome, with center 

marker 


R*B + (S mod B) 
— > R 


S/B—>S 


Sl PT 


The logic of a 
subroutine to 
establish the 
palindromic nature 
of a number. 


--D.J.B. 1979 


247 187 135 91 375 374 373 372 298 400 370 369 368 367 112 160 216 280 352 
296 188 136 92 302 301 300 299 233 380 325 324 295 294 293 159 215 279 232 
230 229 137 93 57 236 235 234 176 306 379 326 257 256 228 158 214 175 174 


172 171 170 94 58 179 178 177 127 240 305 378 327 258 197 196 126 125 124 
123 122 243 212 276 59 31 129 128 86 182 239 304 377 328 259 85 84 83 82 
386 511 244 275 347 395 32 88 87 53 132 181 238 303 376 52 51 56 285 356 
387 312 274 346 396 320 252 13 54 28 90 131 180 237 27 26 166 221 284 355 
388 313 345 397 321 253 193 141 29 11 56 89130 10 79 118 165 220 283 354 
389 344 398 322 254 194 142 98 62 2 30 55 24 47 78 117 164 219 282 353 
390 399 323 255 195 143 99 63 35 15 12 6 8 23 46 77 116 163 218 281 
315 248 189 138 95 60 33 14% 3 1 9 25 49 81 121 169 225 289 361 381 
391 316 249 190 139 96 61 34 4 20 7 48 80 120 168 224 288 360 382 351 
392 317 250 191 140 97 16 5 71 42 21 19 119 167 223 287 359 383 307 350 
393 318 251 192 36 17 18 154 109 72 43 4O 41 222 286 358 384 308 278 349 
394 319 64 37 38 39 269 208155110 73 69 70 22 357 385 309 241 277 348 
100 65 66 67 68 366 339 270 209 156 111 106 107 108 45 310 242 213 144 lol 
102 103 104 105 227 292 365 340 271 210 157 151 152 153 44 76 183 145 146 147 
148 149 150 184 198 226 291 364 341 272 211 204 205 206 207 75 115 199 200 201 
202 203 245 185 133 260 261 290 363 342 273 265 266 267 268 74 114 162 262 263 
264 314 246 186 134 329 330 331 332 362 343 334 335 336 337 338 113 161 217 333 


The end result of the transformations called for 
in the 20 x 20 Problem described on page 2. 

As a coding exercise, the end result would be 
checked sufficiently by printing just the main diagonal 
elements of the array: 

A(i,1) = 371 
A(2,2) = 296 


A(20,20) = 333 


